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1.  Introduction  v 

The  present  work  is  concerned  with  the  determination  of  the  effec¬ 
tive  elastic  properties  of  a  body  containing  many  interacting  flat 
cracks.  The  problem  is  of  significant  practical  importance.  For 
example:  it  has  been  observed  that  crack  density  in  rocks  increases 

before  the  occurrence  of  an  earthquake.  If  it  were  possible  to  determine 
the  stiffness  reduction  due  to  cracks,  analytically,  the  measurement  of 
rock  elastic  moduli  could  provide  information  about  the  density  of  crack 
accumulation.  Another  example  concerns  fatigue  damage.  In  unidirectional 
fiber  composites  fatigue  damage  frequently  accumulates  in  the  form  of 
cracks  parallel  to  to  fibers.  This  reduces  the  elastic  moduli  of  the 
composite,  a  process  which  has  become  known  as  "wearout".  The  crack 
density,  thus  a  measure  of  the  internal  fatigue  damage,  could  be  easily 
determined  by  elastic  moduli  reduction  measurement  if  a  quantitive  re¬ 
lation  between  moduli  and  crack  density  were  available. 

The  problem  under  consideration  has  received  repeated  attention. 

Early  work  was  primarily  concerned  with  the  case  of  low  crack  density 
which  implies  that  cracks  are  sufficiently  far  removed  from  one  another 
so  that  the  stress  field  around  each  can  be  approximated  by  the  stress 
field  around  an  isolated  crack  in  an  infinite  medium.  The  first  work 
of  this  nature  appears  to  be  that  of  Bristow  [1]  who  treated  the  plane 
case  of  long  rectangular  aligned  cracks  extending  indefinitely  on  x^ 
direction  and  randomly  oriented  in  the  Xj,  x2  plane.  Walsh  [2,3,4]  solved 
the  case  of  randomly  oriented  plane  elliptical  cracks,  Salganic  [5]  re¬ 
considered  the  same  problem  by  a  dislocation  approach.  Knopoff  and  Garbin 
[6,7]  treated  this  problem  in  terms  of  scattering  of  long  waves  by  the 
cracks.  Piau  [8]  used  such  an  approach  for  the  problem  of  aligned  cir¬ 
cular  cracks.  It  should  be  pointed  out  that  the  scattering  approach  is 
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very  complicated  and  that  the  same  results  are  obtained  in  much  easier  j 

fashion  by  static  methods. 

The  case  of  interacting  cracks  is  very  difficult.  The  only  exact 
result  available  appears  to  be  due  to  Delameter,  Herrmann  and  Barnett 
[9]  who  treated  the  case  of  a  periodic  rectangular  plane  array  of  cracks 
in  terms  of  periodic  dislocation  fields  and  solved  the  resulting  integral 

equations  numerically. 

/ 

The  self  consistent  scheme  (SCS)  approximation  has  been  applied  by 
Budiansky  and  O'Connell  [10]  to  assess  the  effective  elastic  properties 
of  an  isotropic  body  containing  randomly  oriented  elliptical  cracks. 

The  same  method  was  used  by  Hoenig  [11]  who  considered  oriented  circular 
cracks  in  isotropic  material.  Levin  [12]  gave  an  approximate  treatment 
using  the  Green  tensor,  assuming  that  the  influence  of  any  crack  on  the 
other  cracks  can  be  replaced  by  the  influence  of  a  force  dipole  situated 
at  the  crack  center. 

In  view  of  the  great  analytical  difficulty  involved  in  direct  computa¬ 
tion  of  effective  moduli  the  problem  is  here  considered  primarily  by 
variational  methods  which  permit  establishment  of  bounds  on  the  effective 
moduli.  Here  we  shall  also  be  concerned  with  the  case  of  orthotropic 
matrix  which  is  of  importance  for  fiber  composites.  All  of  the  available 
treatments  in  the  literature  seem  to  be  confined  to  isotropic  matrix. 

2 .  Formulation 

An  elastic  body  contains  a  large  number  M  of  volumeless  flat  stable 
cracks  whose  distribution  is  statistically  homogeneous.  The  external 
surface  S  is  subjected  to  the  homogeneous  boundary  conditions  j 
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Ti(s'  =  a'f.n.  =  TV  (2.1) 

whore  are  constant  stresses.  In  the  absence  of  cracks,  is  the 

stress  field  in  the  body.  In  the  presence  of  cracks  the  stress  field 
becomes  very  complicated  since  it  must  in  addition  satisfy  the  conditions 

VV  =  0  (2.2) 

on  all  crack  surfaces  S  .  However,  by  the  average  stress  theorem,  [13], 


°ij  =  °ij  (2-3) 

where  here  and  from  now  on  overbar  denotes  average  over  the  entire  body, 

and  by  statistical  homogeneity  also  over  any  representative  volume  ele¬ 
ment  (RVE)  containing  sufficiently  many  cracks. 

The  effective  elastic  properties  of  a  composite  material  or  hetero¬ 


geneous  medium  are  defined  by  the  relations 

* 


a  ■  ■  =  C  •  - ,  ,  e 
ii  ilkl  kl 

(a) 

* 

eij  =  Sijkl  Gkl 

(b) 

(2.4) 


where  and  are  effective  moduli  and  compliances,  respectively, 

arc  given  by  (2.3)  and  are  the  average  strains.  It  should  be 
noted  that  for  a  material  containing  holes  or  cracks  are  average 

strains  taken  over  matrix  and  holes  as  if  the  latter  were  deforming 
continua. 

Decomposing  average  strains  and  stresses  into  weighted  averages  in 
terms  of  phase  region  averages  and  their  volume  fractions  and  using  the 
elastic  matrix  stress  -  strain  relations  it  follows  that 

S*  =  S(1)  d  +  ?  (2)v 

l j  kl  b  ijkl  °kl  cij  v2 


(2.5) 
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where  1  denotes  matrix  and  2  denotes  holes  of  any  shape.  The  average 

-  ( 2 1  • 

strain  theorem  (see  e.g.  [13])  can  be  used  to  transform  e  —  into  a 
surface  integral  over  the  holes.  In  the  special  case  when  the  holes 
become  flat  volumeless  cracks  (2.5)  assumes  the  form 


Jijkl  °k] 


+  Yii 


(a) 


(2.6) 


Y.  . 

l  I 


?V 


7  l  f  (I 

s 

m 


u .  In  , 


m 


+  f  u  .  1  n  . )  tlS 
.1  i 


(b) 


where  Sm  is  the  surface  of  the  mth  crack  and  [u^]  are  the  displacement 
jumps  across  the  crack  surfaces. 

If  the  body  is  statistically  isotropic  with  isotropic  matrix  it 

* 

follows  from  (2.6)  that  the  effective  bulk  modulus  K  and  the  effective 

* 

shear  modulus  G  are  given  by 
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(2.7) 


where  o°  is  an  applied  average  isotropic  stress  and  is  an  applied 

a/erage  shear  stress.  Expressions  of  type  (2.7)  can  easily  be  written 
down  for  cases  of  anisotropic  matrix  and  statistical  anisotropy  dic¬ 
tated  by  crack  distribution. 

If  follows  that  the  effective  elastic  properties  of  cracked  bodies 
are  determined  by  the  crack  opening  displacements.  To  determine  these 
it  is  in  general  necessary  to  find  the  entire  displacement  field  in  the 
cracked  body  subject  to  (2.1)  and  (2.2).  This,  of  course,  is  an  enor¬ 
mously  difficult  problem. 
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Effcctive  elastic  moduli  of  heterogeneous  materials  can  alterna 
tively  and  equivalently  be  defined  in  terms  of  elastic  energy.  The 
stress  energy  stored  in  any  heterogeneous  body  subject  to  (2.1)  is 
rigorously  given  by,  [13]. 
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(2.9) 


and  6Um  is  the  energy  change  due  to  the  m  crack  in  the  presence  of 
all  the  others. 

According  to  a  well  known  result,  Eshelby,  [14],  6Um  for  a  hole 
can  be  written  as 

{ a«  =  ?  Kj  “iV s  (2’10) 
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which  in  the  case  of  a  crack  becomes 
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Eqns.  (2.8-11)  may  be  summarized  in  the  form 
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J_ 
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.Y.  V 


(2.12) 


which  is  the  equivalent  of  (2.6). 

Alternatively  (2.11)  can  be  expressed  in  terms  of  stress  intensity 
factor  (SIF)  integrals.  For  line  cracks  in  plane  isotropic  elasticity 


■\i.n  =  —  ju^(x)  +  :<^(x)  +  ( ,_vUni(x)J'JX  (2.13) 

o 

2 

where  =  1  for  plane  stress,  <J>  =  1-  v  for  plane  strain,  am  is  half  the 
crack  length  and  Kj,  Kjj,  K  ^  ^  j  are  SIF  in  modes  I,  II  and  III  respec¬ 
tively.  For  an  isolated  crack  in  an  infinite  body  (2.13)  assumes  the 
form 
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a  Uf 
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II 


(1-V),<I1I  1 


(2.14) 


Corresponding  results  ace  available  for  plane  cracks  of  arbitrary  plane 
shapes  and  in  particular  for  elliptical  cracks,  [10]. 

5 .  Direct  Approaches 

The  case  which  is  most  easily  solved  is  low  crack  density.  In  this 
event  it  is  assumed  that  cracks  are  sufficiently  far  apart  so  that  the 
energy  change  due  to  any  one  crack  is  determined  as  if  it  were  isolated 
in  infinite  matrix  with  the  applied  homogeneous  stress  field  at  infinity. 

Examination  of  the  form  of  5U  for  plane  and  elliptical  cracks  reveals 
that  for  low  crack  density  (2.8-9)  are  equivalent  to 
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where  depends  on  matrix  properties  and  crack  geometry  and  the 

crack  density  parameter  is  given  by 
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(3.2) 


where  in  the  first  case  am  is  half  the  crack  length  and  A  is  the  area 
of  the  plane  body  while  in  the  second  case  a  ,  b  are  major  and  minor 
axes  of  the  elliptical  crack  respectively,  while  V  is  the  volume.  Since 
oV  and  0°^  are  arbitrary  stresses  it  follows  that 


°i jkL  '  Jijkl 
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(3.3) 


where  ~  denotes,  here  and  from  now  on,  low  crack  density  result. 

As  an  example  we  consider  a  plane  orthotropic  sheet  containing 
cracks  aligned  in  direction,  in  a  state  of  plane  average  stress  o°^, 
°22’  °12’  The  strain*  stress  relations  of  the  uncracked  material  are 
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This  may  serve  as  a  model  for  a  unidirectionally  fiber  reinforced  lamina 
in  which  case  Ej  =  E^  -  axial  Young’s  modulus;  E2  =  E^,  -  transverse  Young's 
-  axial  Poisson's  ratio;  G12=GA 


modulus;  V^2=VA 


axial  shear  modulus. 


1 ne  change  in  stress  energy  due  to  a  crack  of  length  2a^  in  an  orthotropic 
sheet  is  given  by,  (15,  10] 

AUBsiiTt{  (3.5) 
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For  an  isolated  crack  in  an  infinited  sheet  this  becomes 
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(3.7) 


Since  the  crack  does  not  perturb  a  homogeneous  o^j  field  this 
stress  does  not  enter  into  the  energy  change.  Comparison  of  (3.7)  with 
(3.1)  and  (3.3),  taking  into  account  the  orthotropy  of  cracked  and  un¬ 
cracked  material  yields 
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where  u  is  given  by  (3.2a). 

It  should  be  noted  that  (3,9)  are  exact  results  for  any  aligned 
crack  distributions  since  such  cracks  do  not  perturb  a  homogeneous 
stress . 

When  the  crack  density  is  not  low  and  thus  crack  interaction  cannot 
be  neglected  the  problem  becomes  extremely  difficult.  In  view  of  the 
formulation  given  above  it  is  necessary  to  either  know  the  crack  face 
displacements  or  their  stress  intensity  factors.  To  obtain  such  informa¬ 
tion  it  is  in  general  necessary  to  determine  the  displacement  fields  or 
the  singular  stress  fields  in  the  cracked  body.  The  only  exact  treat¬ 
ment  available  seems  to  be  due  to  Delameter,  Herrmann  and  Barnett  [9] 
for  the  case  of  a  rectangular  periodic  array  of  plane  cracks  in  isotropic 
matrix.  Using  dislocation  theory  the  problem  was  formulated  in  terms 
of  integral  equations  which  were  solved  numerically.  ([9b]  corrects 
numerical  mistakes  appearing  in  [9a]). 

An  approximate  method  known  as  the  self  consistent  scheme  (SCS)  has 
first  been  applied  to  the  problem  by  Budiansky  and  O'Connell  [10].  The 
basic  idea  is  very  simple:  it  is  assumed  that  the  energy  change  associ¬ 
ated  with  any  crack  can  be  estimated  as  if  this  crack  were  located  in  a 
homogeneous  elastic  body  whose  elastic  properties  are  the  effective 
elastic  properties  to  be  determined.  The  analytical  consequence  of  such 
an  assumption  is  as  follows:  if  the  dependence  of  in  (3.3)  on 

matrix  elastic  moduli  C • • , .  or  compliances  S  .  . ,  .  is  denoted  r..,.(C)  or 

ljkl  r  likl  likl  — 1 
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(3.10) 


where  are  the  effective  compliances  to  be  determined.  Equs. 

(3.10)  are  a  set  of  simultaneous  equations  to  determine  the  effective 
compliances . 

While  this  approach  is  simple  and  straightforward  it  must  be  realized 
that  it  is  based  on  the  drastic  assumption  that  the  crack  tip  "sees" 
the  effective  medium.  This  is  of  course  incorrect  since  the  crack  tip 
sees  matrix  and  other  cracks.  The  concept  of  the  effective  medium  is 
valid  only  for  averages  over  RVE  containing  many  cracks.  It  can  certain¬ 
ly  not  be  applied  to  the  immediate  vicinity  of  a  crack  tip  where  the 
stress  gradients  are  extremely  large.  The  energy  change  due  a  a  crack 
may  be  accurately  estimated  by  the  SCS  method  when  the  crack  is  by  an 
order  of  magnitude  larger  than  all  others,  but  not  necessarily  when  all 
are  of  same  order  of  magnitude. 

As  an  example  for  (3.10)  we  consider  again  the  cracked  orthotropic 
sheet  with  cracks  aligned  in  x^  direction.  Combination  of  (3.10)  with 
(3.8)  gives  the  equations 


(3.11) 


Eqs.  (3.11)  are  equivalent  to  $  order  algebraic  equations  for 
each  of  the  effective  moduli  and  therefore  have  multiple  real  roots. 
Solution  of  the  equations  for  matrix  properties  typical  of  unidirectional 
glass  epoxy  fiber  composite  gave  two  different  £2(01)  both  mononotonically 
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decreasing  with  a  and  two  different  G^C01)  one  which  decreases  with 

* 

a  while  the  other  increases  with  a.  The  second  G^C01)  is  therefore 

* 

physically  unacceptable.  To  make  a  choice  between  the  two  F?(ct)  it 

would  be  necessary  to  argue  that  the  one  associated  with  unacceptable 
* 

G^(u)  via  (3.11)  is  also  unacceptable.  Whether  or  not  such  an  argument 
is  convincing  remains  an  open  question. 

It  is  of  interest  to  note  that  similar  problems  arise  with  the  SCS 
solution  for  aligned  circular  cracks  given  by  Hoenig  [11],  but  only  a 
single  numerical  solution  has  been  presented  in  his  paper. 

Finally  it  is  noted  that  the  SCS  may  easily  be  generalized  to  the 
case  when  there  are  aligned  cracks  in  and  directions.  The  results 
which  have  been  given  in  [17]  are: 

^S'V^'V'7'  -  V^/J^Y  -  2v1?)-|<i-n,)™2l1/2  =  o 
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and  a.  and  a,  are  half  crack  lengths  of  cracks  in  x,  and  x?  directions 
lm  2m  1  L 

respectively.  These  equations  also  have  multiple  real  roots. 


4 .  Bounds  for  Effective  lilastic  Moduli 

The  problem  of  bounding  the  effective  elastic  properties  of  cracked 
materials  will  here  be  considered  in  terms  of  the  classical  extremum 
principles  of  minimum  potential  and  complementary  energies. 

Upper  Bounds 

To  construct  upper  bounds  we  shall  use  the  principle  of  minimum 
potential  energy.  The  potential  energy  and  the  potential  energy  func¬ 
tional  are  defined  by 


(a) 

(4.1) 

(b) 

where  W  is  the  strain  energy  density,  T^  are  prescribed  tractions  on  ST, 
ru  is  an  admissible  displacement  field  which  must  be  continuous  in  the 
region  excluding  the  cracks  and  satisfy  displacement  boundary  conditions, 
if  imposed,  and 
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The  extremum  principle  states  that 


U  >  U 
P  -  P 


(4.3) 


In  the  case  of  a  cracked  body  it  follows  from  (2.2)  that  the  crack 
surfaces  do  not  contribute  to  the  surface  integral  in  (4.1).  Thus  S-j. 
becomes  the  external  surface  S  on  which  (2.1)  is  prescribed.  Since 
o°j  are  constant  it  follows  at  once  that 


-  1  (  o  -  . )  c.  d‘J 

p  '  T  ij  ij  iJ 


(4.4) 


\  -■  )f  ‘"u  - 


The  last  equation  (4.4b)  following  from  the  average  theorem  of  virtual 
work,  [13],  and  (2.4b).  It  is  noted  in  passing  that  (4.4)  apply  for 
any  heterogeneous  body  with  (2.1)  prescribed. 

An  admissible  displacement  field  for  a  body  containing  an  arbitrary 
distribution  of  cracks  is  here  chosen  in  the  form 


'V  o  V  , 

U  .  ;  U  4-  l  u! 


i  1 


(4.5) 


whe  re 


o  o  „  o  o 

u=e.  x.  e. 
i  ij  J  ij  ij*!  <1 


(4.6) 


is  the  displacement  field  in  the  absence  of  cracks  and  u’  is  the  pertur- 
r  im  r 

t  h 

bation  field  of  the  in  crack,  as  if  all  the  other  cracks  were  absent. 
Thus  for  cracks  removed  from  the  boundary  u°+u|m  is  the  displacement 
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ficld  when  the  crack  is  isolated  in  an  infinite  medium.  The  field 
(4.5)  is  the  actual  solution  for  low  crack  density.  It  is  continuous 
in  the  region  excluding  the  cracks  and  is  thus  an  acceptable  admissible 
field  for  that  region.  Across  each  crack  (4.5)  is  discontinuous  with 
jump  [u1]  of  the  displacement  of  the  isolated  crack.  It  should  be  noted 
that  an  admissible  displacement  field  with  no  discontinuity  across 
cracks  would  lead  to  the  trivial  result  that  the  upper  bounds  for  ef¬ 
fective  moduli  are  the  matrix  moduli. 

The  strains  and  stresses  associated  with  (4.5)  are 
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Since  o- .  is  the  actual  stress  perturbation  due  to  an  isolated  crack 
l  j  m  1 

with  (2.1)  prescribed,  each  of  the  associated  tractions  vanish  on  the 
external  boundary.  Thus 


r  =  o .  .  n .  =  n 
im  iJ»>  I 
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(4.8) 


Introduction  of  (4.7)  into  (4.1b)  yields  the  result 
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lt  is  easily  shown  that  the  term  in  the  first  integral  of 

(4.9a)  is  -  6Um,  the  negative  of  the  stress  energy  change  due  to 
an  crack  isolated  in  infinite  medium.  Therefore  (4.9)  assumes 
the  form 
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The  stress  o°.  +  o 
satisfies  equilibrium, 
equilibrium.  Similarly 
by  virtual  work 


is  an  actual  stress  field  and  therefore 

ljm 

o  * 

Since  a-.  is  constant,  o..  also  satisfies 
ij  ljm 
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e..  is  a  compatible  strain  field.  Therefore 
i  j  n  r 
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on  all  cracks 

The  first  integral  on  the  right  hand  side  vanishes  because  of  (4.8). 
The  second  integral  must  be  evaluated  on  the  two  adjoining  faces  of 
all  cracks  on  which 
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Thus  the  integrand  becomes  -  o..u.  n-  and  since  u.  is  continuous  across 
6  ij  in  j  in 

all  cracks  and  since  the  signs  of  the  normal  are  opposite  on  adjoining 
crack  faces,  the  last  integral  in  (4.11)  also  vanishes  and  thus  (4.11) 
is  zero.  It  then  follows  from  (4.10)  that 


0P  =  -H  -  l  (S  || 

1  o  L  n 

in 


(4.12) 


Now  6Um  is  the  energy  change  due  to  one  isolated  crack  in  a  large  body. 
Most  cracks  are  sufficiently  removed  from  the  boundary  so  that  6Um 
can  be  computed  as  if  the  crack  were  at  infinite  distance  from  the 
boundary,  thus  by  expressions  of  type  (2.14).  The  cracks  which  are 
near  the  boundary  constitute  the  usual  boundary  layer  of  a  hetero¬ 
geneous  body.  Their  number  relative  to  that  of  cracks  removed  from 
the  boundary  can  be  made  indefinitely  small  by  increase  of  the  size 
of  the  cracked  body.  Therefore  (4.12)  can  be  determined  to  any  de¬ 
sired  accuracy  by  interpreting  all  5Um  as  energy  changes  of  isolated 
cracks  removed  from  the  boundary.  It  should  be  noted  that  the  same 
reasoning  tacitly  enters  into  derivation  of  all  previous  results  such 
as  low  crack  density,  periodic  crack  arrays  and  SCS. 

If  follows  from  the  low  crack  density  analysis  leading  to  (3.1), 
(3.3)  that  (4.12)  can  be  written  as 


ur  =  - 


1  *  o  o 
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Then  from  (4.3)  and  the  last  of  (4.1b) 


(4.13) 
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(4.14) 


This  positive  definite  form  defines  bounds  for  S  in  terms  of  the 

low  crack  density  expressions  as  given  by  (3.3)  in  which,  how¬ 

ever,  a  is  no  longer  small  but  arbitrary.  Specializing  (4.14)  to  the 
case  of  one  normal  zero  stress,  say,  it  follows  that  S 


2222 


is  smaller 
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than  $2222'  Since  t^e  reciprocals  of  these  compliances  are  Young's 
moduli  in  2  direction  we  have 


,.*(+) 
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(4.15) 


Analogously,  for  one  nonzero  shear  stress,  0°^  say>  it  follows  that 
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(4.16) 


Similar  results  are  valid  for  the  other  Young's  and  shear  moduli  for 
any  cracked  material.  Thus  the  effective  moduli  expressions  val id 
fo_r  low  crack  density  are  upper  bounds  for  th£  effective  modul i 
at  any  crack  density. 

If  the  cracked  material  is  statistically  isotropic  (4.14-16)  be¬ 
come  bounds  for  that  case.  It  is  also  easily  shown  that  the  bulk 
* 

modulus  K  is  bounded  by 


*  A  * 

K  <  K  =  .< 


(4.17) 


For  randomly  oriented  circular  cracks  the  results  are. 


where 
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and  am  is  the  radius  of  the  m  crack.  Thus  (4.18)  are  upper  bounds 
*  * 

for  E  and  G  of  a  body  containing  an  arbitrary  number  of  circular 
cracks . 

For  an  orthotropic  sheet  which  contains  an  arbitrary  number  of 

cracks  aligned  in  direction  it  follows  from  (4.15-16)  that  (3.8) 

*  * 

become  upper  bounds  for  the  moduli  E2  and 

Lower  Bounds 

To  construct  lower  bounds  we  shall  use  the  principle  of  minimum 
complementary  energy.  Since  the  imposed  boundary  conditions  (2.1), 

(2.2)  are  traction  boundary  conditions  the  principle  simplifies  to 
that  of  minimum  stress  energy.  The  stress  energy  and  the  stress  energy 
functional  are  defined  by 


a!  r  1  *  o 
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where  is  an  admissible  stress  field  which  must  satisfy  equilibrium 
and  all  traction  boundary  conditions. 


The  principle  of  minimum  stress  energy  states 

0  i  Ua 


(4.20) 


The  major  difficulty  in  construction  of  is  satisfaction  of 

the  crack  boundary  conditions  (2.2).  Consider  a  body  with  an  arbitrary 

distribution  of  cracks,  fig.  1.  The  body  is  subdivided  into  regions 

each  of  which  contains  one  crack.  Let  it  be  assumed  that  the  m1”*1  such 

region  is  subjected  on  its  boundary  Sm  to  (2.1).  The  elastic  stress 

field  in  this  body  satisfying  (2.1)  and  (2.2)  is  denoted  0?^  Then 
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evidently  the  stress  field 
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(4.21) 


is  admissible  for  the  entire  body  since  it  satisfies  equilibrium  and 
all  boundary  conditions. 

It  should  be  noted  that  the  proposed  a- .  is  singular  at  the  crack 

tips  and  the  usual  admissible  fields  considered  in  standard  elasticity 

texts  do  not  seem  to  include  such  cases.  It  may  be  argued  that  cracks 

are  limiting  cases  of  very  flat  holes  with  finite  tip  curvatures  and 

no  stress  singularities.  Then  the  stress  energy  functional  (4.19b) 

must  first  be  evaluated  for  stresses  in  regions  containing  such 

holes  and  the  limit  of  these  integrals  are  taken  when  the  holes  become 

cracks.  These  integrals  are  the  actual  strain  energies  stored  in  these 

regions  subjected  to  (2.1)  and  each  containing  one  crack  and  it  is 

t  h 

known  that  these  energies  remain  finite.  We  can  write  for  the  m 
region 
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and  6Um  can  be  expressed  in  terms  of  the  stress  intensity  factor  of 
the  crack  in  the  m**1  region,  if  known. 

For  line  cracks  and  plane  elasticity  $Um  will  have  the  form  (2.13). 
It  is  emphasized  that  the  SIF  entering  into  the  expressions  depend 
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on  the  shape  of  the  region  and 

If  follows  that 


its  crack  geometry. 
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A  crack  solution  useful  for  present  purposes  is  due  to  Isida 
[18]  who  considered  the  plane  case  of  a  finite  rectangle  with  a  central 
crack  in  Mode  I.  fig.  2.  His  expression  for  the  SIF  may  be  written 
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where  oQa  is  the  mode  I  SIF  for  the  same  crack  in  plane  stress  when 
it  is  isolated  in  an  infinite  sheet.  Numerical  values  of  the  function 
fj  have  been  given  in  [18]. 

The  same  problem  for  an  orthotropic  rectangle  has  been  treated 
numerically  by  Bowie  and  Freese  [19].  The  SIF  for  this  case  may  be 
expressed  as 

1/2 

=  oQa  f^fa/b.a/c,  material  properties)  (4.25) 

Consider  now  a  large  plane  specimen  containing  m  randomly  dis¬ 
tributed  cracks,  all  in  x^  direction,  fig.  3.  We  can  subdivide  this 
specimen  into  rectangles,  each  rectangle  containing  a  central  crack 
or  no  crack  as  shown  by  the  dashed  lines.  The  aggregate  of  the  solutions 
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of  finite  rectangles  with  dimensions  2b  and  2c  and  central  cracks 

°  mm 

of  length  2a  (a  can  also  vanish)  under  simple  tension  a  is  an 
°  m  m  r  o 

admissible  stress  fields.  In  the  special  case  when  the  ratios  a  /b  , 

r  mm’ 

a  /c  and  b  /c  remain  constant  throughout  the  specimen,  eqs.  (4.21  - 
mm  mm 

* 

24)  together  with  (2.13)  give  a  lower  bound  on  E2  for  isotropic 
matrix  with  aligned  cracks  in  the  form 
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[1  +  f  fj2(a/b,a,c)a  dal 

a  ' 
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(4.26) 


For  orthotropic  matrix  we  use  (4.25)  together  with  (3.5)  and  obtain 

k 

a  lower  bound  on  E  2  in  the  form 

* 

K2 
K 


The  bounds  (4.26)  and  (4.27)  are  evidently  valid  also  for  the  case  of 
a  doubly  periodic  array  of  cracks,  since  in  this  case  a^,  bm  and  cm 
are  constant. 

Another  useful  crack  solution  is  due  to  Wilson  [20],  who  consider 
ed  the  plane  case  of  a  finite  isotropic  rectangle  with  a  central  in¬ 
clined  crack  under  various  loadings. 
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The  solution  for  a  finite  rectangle,  in  tension  or  in  shear, 

containing  an  inclined  crack  can  be  used  to  obtain  lower  bounds  on 

*  * 

E  and  G  for  a  plane  sheet  with  randomly  distributed  cracks.  The 
solution  for  a  finite  rectangle  with  a  central  oblique  crack  sub¬ 
jected  to  uniform  shear  tractions  on  crack  edges  can  be  used  to  ob- 

* 

tain  a  lower  bound  on  G^  for  a  plane  isotropic  sheet  with  aligned 
cracks . 

Some  of  the  results  obtained  are  plotted  in  fig.  4  as  a  function 
of  the  crack  density  parameter  a.  This  parameter  is  the  only  geo 
metrical  information  entering  into  the  upper  bounds  (small  crack  density) 
and  into  the  SCS.  However,  a  is  insufficient  information  to  charac¬ 
terize  the  exact  results  given  in  [9]  or  the  lower  bounds.  In  both 
cases  the  results  are  based  on  specific  geometries  which  are  charac- 

7 

terized  by  the  ratios  c/b  and  a/b.  Then  a.  ~  a  /4bc  and  is  clearly 
insufficient  to  specify  the  geometry.  The  c/b  ratios  used  for  lower 
bound  calculations  are  indicated  in  the  figures. 

The  upper  and  lower  bounds  are  quite  close  for  a  respectable  range 
of  x ,  up  to  a  -  0.15.  For  small  a(a<0.05)  all  the  results  coincide. 

The  SCS  curves  are  close  to  the  lower  bounds.  As  the  cracks  become 
closer,  i.e.  for  large  a/b  and  small  c/b,  the  distance  between  the 
bounds  increases,  and  the  small  concentration  result  deviates  further. 
This  is  reasonable,  as  the  influence  of  the  interactions  between  the 
cracks  should  become  more  important.  The  results  are  presented  for 
a  range  of  up  to  a  =  0.3,  however  the  cracks  may  become  unstable  be¬ 
fore  a  attains  this  value. 

It  has  been  found  that  in  certain  cases  the  results  of  [9]  are 
above  the  general  upper  bound,  which  is  unacceptable.  The  reason  for 
this  phenomenon  is  not  known  to  us  at  the  present  time. 
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5. 


The  problem  of  analytical  determination  of  elastic  properties  of 
cracked  materials  has  been  considered  by  variational  methods.  The 
motivation  for  this  approach  is  that  exact  direct  analysis  of  the 
problem  is  extremely  difficult.  The  only  exact  solution  available  for 
arbitrary  crack  density  is  for  a  periodic  crack  array  and  had  to  be 
performed  numerically. 

The  bounds  obtained  are  reasonable  close  for  crack  densities  of 
practical  significance.  When  the  bounds  are  close  the  small  crack 
density  results  can  be  regarded  as  a  good  approximation  since  it  has 
been  shown  that  they  are  general  upper  bounds. 

On  the  basis  of  the  results  obtained  it  appears  that  the  SCS  ap¬ 
proximation  tends  to  underestimate  the  effective  elastic  properties. 

A  more  serious  problem  with  this  method  is  that  it  can  give  non-unique 
results  for  effective  elastic  moduli. 

The  bounds  derived  in  this  work  are  based  on  known  crack  field 
solutions  and  require  only  small  computational  effort.  It  is  to  be 
expected  that  additional  and  perhaps  better  bounds  can  be  derived  as 
more  crack  solutions  become  available. 
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Fig.  1  -  Subdivision  of  cracked  body  for  construction 
of  admissible  stress  fields:  General  case. 
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ig.  2  -  Finite  rectangle  with  a  central  crack. 
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Bounds  and  results  for  Young's  and  shear  moduli  (a)  Isotropic  matrix 
randomly  distributed  cracks  (b)  Isotropic  matrix,  aligned  cracks 
(c)  Orthotropic  matrix  (glass/epoxy  composite  =  26.4  GPa; 

E?  =  7.17  GPa;  G17  =  5.12  GPa;  v12  =  .267)  with  aligned  cracks. 


